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ABSTRACT 



I describe a unified formalism for lattice fermions, in which the relativistic 
action of Wilson and the nonrelativistic and static actions appear as special 



(N ' 

on : 

Q\ , cases. It is valid at all values of m q a, including m q a ~ 1. In the limit 

m q a <C 1, the formulation reduces to the light quark action of Wilson. In 
^ ■ the limit m q a 3> 1, the formulation reduces to the nonrelativistic action of 

Q ■ Thacker and Lepage, and to the static action of Eichten. 

-xt- : 

Present and future lattice calculations involving b and c quarks include some of the 
CNJ ■ most important applications of lattice gauge theory to standard model physics. These 

include the heavy meson decay constants, the BB mixing amplitude, and various 
,_i , semileptonic decay amplitudes, which are all crucial in extracting CKM angles from 

experimental data. They also include the extraction of a s from the charmonium and 
| bottomonium spectra. 

CN ■ There exist two main classes of methods for lattice fermions: small mass ex- 

pansions and large mass expansions. The standard Wilson action contains errors 
which vanish as powers of m q a, the quark mass in lattice units. This series of cor- 
rections does not converge when m q a > 1. There exists a second class of methods 
for treating lattice fermions when m q a > 1: Nonrelativistic QCD (NRQCD)S and 
the static approximation.! Because coefficients of higher terms in the Lagrangians 
of these methods (such as D 2 /(2m)) are explicit functions of 1/m, the loop correc- 
^ tions are also explicit functions of 1/m. These begin to diverge as ma is reduced 

below a value of order one, making the nonrelativistic expansion impractical. The 
masses of the b and c quarks are such that m q a is 0(1) at the lattice spacings used 
in current numerical work. Therefore, calculations of such crucially interesting quan- 
tities as the heavy meson decay constants and fr> have often involved awkward 
interpolations between results in the static approximation and results using Wilson 
fermions through a region where neither approximation is well behaved. While such 
an approach is probably workable, it is clearly desirable to have a method for lattice 
fermions which is well behaved throughout the region of interest. 

Since the large mass and small mass formalisms are both descriptions of QCD, it 
is not surprising that they share certain fundamental features. Like Wilson fermions 
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(if)), the fermions of NRQCD contain four components per site: a two-component 
quark field (0) and a two-component antiquark field (x). Further, the two actions 
employ the same sorts of interactions: covariant time derivatives, covariant Lapla- 
cians, etc. To find a formalism uniting the two actions, we therefore consider the 
following generalized Lagrangian:i | 

£ = (f)*( ci A-+mo-|E A , +A D0n + c 3 0*E^ A ^n 

1 i i 

+ X *(-ci A+ + m Q -|^A J +A7)xn + c 3 x*E^ A A- (1) 

i i 

(A + , A - , and A are the forward, backward, and symmetric discrete difference opera- 
tors, respectively. I will use m q for the physical quark mass and mo for the bare quark 
mass on the lattice.) With the choice of parameters C\ = c 2 = C3 = 1, this is simply 
the standard Wilson action. When c\ — 1 (times a correction factor when ma 3> 1), 
C2 = — , and C3 is negligible, it is a good, if somewhat unconventional, Lagrangian for 
NRQCD. The bare mass is conventionally omitted in NRQCD calculations, but we 
are free to leave it in the theory. The usual Dirac coupling between quarks and anti- 
quarks is absent (having been transformed into higher derivative interactions by the 
Foldy-Wouthuysen transformation), but we may add back a sufficiently suppressed 
amount of this interaction without spoiling the nonrelativistic theory. It is thus pos- 
sible to adjust the parameters for this particular NRQCD action in such a way that 
as mo is reduced, instead of blowing up, the theory turns smoothly into the Wilson 
theory. One must find appropriate normalization conditions to determine the param- 
eters of the action which lead to the Wilson theory for am q <C 1, and to NRQCD and 
the static approximation for am q <C 1, and which work at all values of am q . 

It is illuminating in this regard to expand the equation for Wilson propagators 
nonrelativistically when the mass is large, to see what breaks down. After normalizing 
the fields by ^==^ 1 (not ^7= as is commonly used) one may obtain 







-S + M + (l- C4, ) - ~ ( — + 1 ) E(a, 




<f>n, (2) 



where S is the energy eigenvalue obtained from the transfer matrix and M. = S p 2 =0 = 
ln(l + mo). This is a lattice Schrodinger equation not unlike the one obtained from 
NRQCD, but it has some unusual features. Most important, the two "masses" in 
the equation, the rest mass M. = ln(l + mo), and the mass governing the energy- 
momentum relation -L = ^ + ( 1+mo ) 1 ( 2+mo ) > are completely different. M. plays 
little dynamical role in heavy quark systems and is usually omitted from NRQCD 
and static approximation calculations. The dynamically more important condition 
dS/dp 2 = l/(2m q ) is used to fix the mass in NRQCD. This condition is also the same 
as the usual mass condition for Wilson fermions when am <C 1, since when m a is 
small, 

1 1 11 
I ~ ~ (3) 

m (1 + m ) (2 + m ) M m 



However, if the rest mass M. is used to fix the fermion mass for Wilson fermions when 
am > 1, the energy- momentum mass condition dS/dp 2 = l/(2m q ), which is usually 
more important, will be completely incorrect. The two masses can be put back into 
agreement with the use of the Lagrangianll I 

C = -Tp n ip n + «tV>„(l - 7o)£4,oVWo + h - c - 

+ k s Y^ $n(l - ^i)U n ^ n+i + h.c. (4) 

i 

Thus, it seems that an action closely related to the Wilson action is a member 
of the class of actions suitable for NRQCD. Further, in NRQCD and in the static 
approximation, M. plays little dynamical role. It can usually be ignored, and is 
usually omitted. This suggests that for that majority of calculations in which the 
mass gap between states containing different numbers of quarks is unimportant, the 
standard Wilson action itself can be used when am > 1 as long as Ai is ignored and 
d£ /dp 2 = is used to fix the quark mass, as is done in NRQCD. 

This proposal is obviously correct in free field theory, where we can calculate the 
behavior of quark propagators exactly to see that the proposed interpretation makes 
sense. It is certainly correct in mean field theory, too. Mean field improvement of 
these fermions, as of Wilson fermions, is simply the absorption of a "mean link" uq 
into an effective k = u k and then proceeding as with free field theory.! (A plausible 
estimate of the mean link in this context is probably uq ~ 1/8k c .) It remains to be 
shown whether the theory is somehow spoiled by renormalization. 

Perturbatively, Green functions must be expanded in p 2 and a s . The coefficient 
of each term in the expansion is an explicit function of the quark mass, since the 
theory must be solved exactly in ma. (The is also the case for the loop corrections 
of NRQCD.0) If these functions become singular or badly behaved in some way, the 
theory could conceivably break down. The one loop perturbative corrections contain 
all of the ugliest features of Wilson and NRQCD perturbation theory simultaneously, 
and have only been begun. There is, however, one numerical calculation by El- 
Khadrai indicating that nothing too surprising occurs. The one-loop correction to 
the local current normalization for Wilson fermions with the naive normalization isi 

= 2.(1 -0.173*) ■ (5 > 
The correct normalization with mean field improvement is 



The remaining perturbative correction, 0.06g 2 , becomes an explicit (so far uncalcu- 
lated) function of m (or k) in the new formalism. The small mass limit of this function 
is 0.06g 2 . This function must not become singular if the theory is to make sense. This 



normalization is straightforward to calculate numerically. As described in Ref. I, the 
nonperturbative calculation agrees with Eq. ^| to within a few per cent. It disagrees 
with Eq. |5| by around a factor of two. It can be seen that for this quantity, not only 
is the unknown function of m not singular, it is approximately equal to 1. 

Putting the new action on a secure footing will ultimately require: 1) determina- 
tion of the bare parameters of the action with mean field theory and full perturbation 
theory, 2) nonperturbative tests of the perturbative results, and 3) phenomenological 
tests of the resulting action in calculations of well understood physical quantities. Not 
much of this program has yet been accomplished. However, from what is known now, 
there appears to be no insuperable obstacle to developing a formulation of lattice 
fermions practical for all values of the quark mass, not just m q a ^> 1 or m q a <C 1. 

Care will clearly be required in formulating normalization conditions which cap- 
ture the most important physics in both the relativistic and nonrelativistic regions. 
(Identifying dS/dp 2 rather than M. as the fundamental mass condition is example 
number one of these.) 

I thank G.P. Lepage and my collaborators A. X. El-Khadra and A. S. Kronfeld 
for useful discussions. Fermilab is operated by Universities Research Association Inc. 
under contract with the U.S. Department of Energy. 

This material appeared in slightly altered form in Ref. |. 
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